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Non-Fermi liquid behavior of the drag and diffusion coefficients in QED plasma 
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We calculate the drag and diffusion coefficients in low temperature QED plasma and go beyond 
the leading order approximation. The non-Fermi-liquid behavior of these coefficients are clearly 
revealed. We observe that the subleading contributions due to the exchange of soft transverse 
photon in both cases are larger than the leading order terms coming from the longitudinal sector. 
The results are presented in closed form at zero and low temperature. 


I. INTRODUCTION 


It has been known for quite sometime now that a fermionic system interacting via the exchange of transverse 
gauge bosons exhibit deviations from the normal Fermi liquid behavior. Such a characteristic feature, in presence 
of transverse or magnetic interactions, for the first time was reported in where the specific heat of a degenerate 
electron gas was shown to contain correction terms involving a,7InT~!. This was interpreted to be a consequence 
of the long range behavior of the magnetic interaction due to the absence of magnetostatic screening. Initially, such 
corrections were considered to be of little practical importance as, such tiny an effect, was not likely to be detected 
experimentally. A decade later, however, the scenario changed and such investigations started attracting attentions in 
the context of strongly correlated electron system in which the gauge coupling is not the fine structure constant (1/137) 
but of order unity |2]. Further impetus to these studies now come from another domain involving relativistic quark 
(or electron) gas at high density and zero or low temperature where the specific heat also contain such anomalous 
terms. This seems to have serious implications in determining the thermodynamic and transport properties of the 
quark component of neutron or proto neutron stars, viz. entropy, pressure, specific heat, viscosity etc . Similar 
non-Fermi-liquid terms, for ungapped quark matter, also appear in the calculation of neutrino emissivity and its mean 
free path (7, . For quarks in a color superconducting state also, the chromomagnetic interaction strongly influences 
the magnitude of the gap as shown in Eg. 

It is well known that the magnetic interaction in non-relativistic systems is suppressed in powers of (v/c)?. The 
scenario, however, changes as one enters into the relativistic domain, where it becomes important. Hence, in dealing 
with relativistic plasma one has to retain both electric and magnetic interactions mediated by the exchange of longi- 
tudinal and transverse gauge bosons like photons or gluons. More interestingly, it is observed that for ultradegenerate 
case, both in QCD and QED, the transverse interactions not only become important but it dominates over its lon- 
gitudinal counterpart; a characteristic behavior having a non-trivial origin residing in the analytical structure of the 
Fermion-self energy close to the Fermi surface. This has been beautifully exposed in where the authors calculate 
Fermionic dispersion relations in ultradegenerate relativistic plasmas and show how such non-Fermi liquid behavior 
emerges from the vanishing of the Fermion propagator near the Fermi surface by calculating the group velocity of 
the corresponding quasiparticle excitations. One can also see, how the fractional power appears there mt due to the 
exchange of soft transverse gauge boson in the small temperature expansion of the Fermion self-energy in ultradegen- 
erate plasma similar to what one encounters in the expansion of the thermodynamic potential or C, (SHB). Actually, 
the fermion self-energy close to the Fermi surface receives a logarithmic enhancement due to the exchange of magnetic 
gluons ad, this, in turn, leads to such dominance. A more rigorous discussion on how and why the dynamics changes 
near the Fermi surface leading to the break down of Fermi liquid behavior or vanishing of the step discontinuity can 
be found in u3. Departure from the Fermi liquid behavior has also been witnessed in the calculation of quasiparticle 
damping rate in ultradegenerate relativistic plasma 

The low temperature, high density region, commonly known as ultradegenerate plasma, is much less explored in 
comparison with the high temperature low density domain. In particular in the present work, we calculate fermionic 
drag (7) and longitudinal diffusion coefficients (B\ ) in this regime and eventually extend it to the limiting case of 
zero temperature. The salient feature here has been the inclusion of the higher order terms both in the transverse 
and longitudinal sector with implications to be discussed later. The evaluation of drag (diffusion) coefficient is very 
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similar to the damping rate calculation with one difference i.e. here we weight the imaginary part of the self energy 
with the energy (square momentum) transfer per scattering to obtain the desired result. Such calculations, as is 
well known, are plagued with infrared divergences. There are well established techniques to handle such divergences 
both at finite and zero temperature where one divides the interactions into two regions one involving the exchange 
of soft gauge bosons while the other involves hard momentum transfer (17). For the former one uses the bare photon 
(gluon) propagator and for the latter the hard thermal/density loop (HTL/HDL) resummed propagator is used. One 
interesting departure from the high temperature that is observed in dealing with plasma close to zero temperature is 
the following: in a hot plasma both the hard and the soft part of the electric and magnetic interactions contribute 
at same order of the coupling parameter. In the ultradegenerate plasma, or when the temperature is much smaller 
compared to the chemical potential, it is seen that the hard sector contribution come with higher order coupling 
parameters than the soft sector. Even within the soft sector, for the longitudinal and transverse part, the coupling 
parameter appears with different powers ig. 

The drag coefficient, as we know, is related to the energy loss suffered by the propagating particle in a plasma. This 
has been studied extensively in a series of works for the last two decades of). There also exist many calculations for 
the diffusion coefficients both for quantum electro and chromomagnetic plasma {19}, p3, B4, Bg). All these calculations 
are performed in situations where the temperature is high but the chemical potential is zero except in Rg, where 
numerical estimates of the energy loss or drag and diffusion coefficients at non-zero chemical potential have been 
presented. There, to the best of our knowledge, exists only one calculation so far {18}, where the analytical results for 
7 and (By ) for ultradegenerate relativistic plasma have been presented. There, we have restricted ourselves only to 
the leading order results and have shown that the drag and diffusion coefficients are dominated by the soft transverse 
photon exchanges while the longitudinal terms are subleading. Here, we go beyond the leading order and reveal 
the importance of the subleading terms in the transverse sector. The approach we adopt in this work is, however, 
different from the previous one and more in line with E. The connections, nevertheless, are made at appropriate 
places. Here, we probably should mention that the dominance of next to leading order (NLO) transverse term over 
the longitudinal one does not imply breakdown of the perturbation series. Because, the next to leading order terms 
in transverse or longitudinal sector individually are smaller than the corresponding leading parts. 

The plan of this paper is as follows: in section II the formalism is set forth. In subsection A, we evaluate the 
drag coefficient in the domain low temperature and eventually arrive at the zero temperature results by taking the 
appropriate limit. In the next subsection we present the results for the diffusion coefficient both at zero and small 
temperature. In section III we conclude. 


II. FORMALISM 
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FIG. 1: Fermion-fermion scattering with screened interaction. 


The drag coefficient of a quasiparticle having energy (E) is incidentally related to the energy loss of the propagating 
particle which undergoes collisions with the constituents of the plasma viz. the electrons: 


n= F | te. (1) 


dI is the differential interaction rate Bi. This expression is quite general and can be used to calculate collisional 
energy loss both for the finite temperature and/or density. The phase space will be different due to the modifications 


of the distribution functions depending upon the values of u and T. The imaginary part of fermion self energy diagram 
basically gives the the damping rate of a hard fermion. This damping mechanism is equivalent to elastic scattering 
off the thermal electrons via the exchange of a collective photon, 


P(E) = -yp [mE (yp +i .p)P|| (2) 


The full fermion self-energy represented in Fig.(I) can be written explicitly as: 
ap) = er | oe one Ss (in — ws), P — A) Ww Apr (iws, q) , (3) 


where, po = iwn + H, qo = iws. Wn = T(2n + 1)T and ws = 2rsT are the Matsubara frequencies for fermion and 
boson respectively with integers n and s. After performing the sum over Matsubara frequency in Eq.(3), iwn + u is 
analytically continued to the Minkowski space iw, + u + po + ie, with € + 0. The blob in the wavy line of Fig. (T) 
represents HTL/HDL corrected photon propagator which is in the Coulomb gauge is given by BI, BÀ, 


Am (Q) = poro Ai(Q) + Pir A:(Q) , (4) 
with, P} = (6:; — âiĝ;), & = at /lal, Ph = Ph; = Pho = 0 and Aj, A, are given by [31 BA, 
=i 
A = — 5 
(Go, 9) +I’ ( ) 
—1 

A = — c. 

(qo, q) q? = g = I; (6) 


Here we introduce the spectral functions pj, as[B1l BA: 


pi t(do; d) 
27 
The poles w; are the solutions of the dispersion relations. The ô function corresponds to the (time-like) poles of the 
resumed propagator and () represent cuts. The latter terms i.e Landau damping pieces of the spectral functions are 
non-vanishing only for qf < q? and are given by, 


mh x O(1 — x”) 


= Zai [6(do — 41,4(@)) — (Go + &1,4(9))] + bildo, a) - (7) 








Bi(qo, 9) =, ioe. S 
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where x = qo/q. The Debye mass is m% = She + = i 


At the leading order these are derived ett the one- fies Siete self-energy where the loop momenta are assumed to 
be hard in comparison to the photon momentum 311, BA. In the literature the formalism is known as the HTL/HDL 
approximation as discussed in [15 

In Eq.(3), fermion propagator has = Dilswine spectral representation with the notation k = (p — q) i 


i dki pr(K 
Sy(iun. 4) = | —— 


(9) 


Taking the imaginary part of Eq. (B), the scattering rate with the help of Eq.(2) can be calculated. One then inserts 
the energy exchange w in the expression of IT and calculate 7 from Eq.(i), to obtain, 


-FJB dko wf E 
1 E j] rP On PEO! J o on 
x (1+ n(qo) — A(ko))6(E — ko — qo) 


[poko + p- k]pi(qo, 4) 
2[poko — (p - â)(k - 4)]px (40,4). (10) 
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The energy conserving delta function in the last equation deletes the contribution from the delta function and therefore 
7 recieves contribution only from the cuts. The same holds true for diffusion coefficients By |; as well. In the above 
equation n and 7 are the Bose-Einstein and the Fermi-Dirac distribution functions: 


1 _ 1 

n(qo) = Ba 1” n(ko) = bto G] (11) 
where, 8 = +. Eq. (5) is the general expression of drag coefficient. 

Apart from 7, the vee baad diffusion coefficient (B;;), could be of i in the study of fermion 


propagating in the plasma p3, p4, Bg. It can be defined as follows [19), p3, 24, 


Decomposing B;j into longitudinal (By ) and transverse components (B1) we get the following expression, 


PiPj PiPj 
Biz = By (di; — mo ) + Bi rE . (13) 








The imaginary part of the Eq. (8) multiplied by the square of the longitudinal momentum transfer in the fermion 
fermion scattering gives the expression for By . Using Eqs. (2) and (13), longitudinal momentum diffusion coefficient 
(Bı = B) can be written as follows, 


Te? d?q f” dko œ dag 
= —— a ——— k eu 2 
B F Qr)? J. on psl of aji 


x (1+ n(qo) — n(ko))d(E — ko — qo) 
x [poko + p - k]pi(qo, @) 
+ 2[poko — (p - 4)(k - 4)]p+(q0, 9). (14) 


Here, q = qcos@ i.e the longitudinal momentum transfer. 


A. Drag coefficient when |E — u| ~ T 


In this section we calculate the drag coefficient (n) when T ~ |E — u| < eu < n, this is the region which is relevant 
for the astrophysical applications. It has been mentioned already that evaluation of 7 is plagued with infrared 
divergences. To circumvent this problem, as mentioned in the introduction, the region of integration as it appears 
below has to be divided into two regions distinguished by the scale of the momentum transfer i.e. the soft and the 
hard sector. For the former, we use the one loop resummed propagator with a finite upper limit on the momentum 
which is designated as q* and for the latter we use the bare photon propagator. Following this prescription, for the 
soft part, one writes: 








8r2E 
x {pi(qx,q) + (1 — x) pr (qx, q)} 


soft e2 q* 3 1 
nf E = Sef dag? f dra + nlar) = a(E ~~ a2) 
0 —1 








(15) 
From the expression after subtracting the energy independent part we have H, 
soft soft e2 q“ 1 
of E-a, = -g |, daa! | dee 0E - u- ar) — van) 
x [(1 — 2) pi(qa, q) + pilaz, q)] - (16) 


First, we calculate the transverse photon contribution then the longitudinal one. For this in Eq. (16) we substitute 
q and qo by introducing dimensionles variables z and v, 


q=2q,2/(nv)'*, qo = Tv, (17) 
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where, qs is the screening distance in the magnetic sector, and we take a = = <1. From the above substitutions it 
immediately follows that, 


2/3 


q= mpa! z, zx =a? vz. (18) 


? 


After expanding the integrand with respect to a we find for the transverse contribution of n, 
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(19) 
where, a = Ezel ~ O(1). Here, we neglect the terms which are more than a™ and G )*. After z integration we 
obtain, 
' soft 25 soft = e?m?, oo i ex —] va? 7 Musas E 20 x 23usas +. (20) 
t t,E=u E Jig (1+e)(1+e%7) \ 24r Ons ITTF 
Now, we use the formula for v integration sending the integration limits to too, 
et — 1 
dy——___—__—_|p|* 
eared 
=T(A +1) [Linga(—e7%) — Ling (-e®)]. VA > 0 
(21) 
soft 
Clearly the expression for 7 is Polylogarithmic in nature, 
t 
soft soft emh a2 
— = — [T (2) (Lie(—e~%) + Lig(—e® 
nl A ea, = EZ) Sag EO) (li-e) + Lia(—e)] 
91/398/3 7 
- 9775 [T($) (Lig/3(—e~%) + Lig/3(—e*)) | 
20 x 2?/3q10/3 ; 2 . a 
= ~o T£) (Liio/3(—€ a + Lijo/3(—e ))] AEren Au (22) 


The above expression can be written in the following form, 
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where, 











From the expression (22) it is evident that the expression contains fractional powers in (E — u). This nature is 
basically a non-Fermi-liquid behavior of ultradegenerate relativistic plasma. After the magnetic part we derive the 
expression of the electric part. 

In case of the electric part we substitute q = qsy and qo = Tu/y, or q = mpy and x = au/y. Though the 
substitutions in electric and magnetic sectors look different, but the nature of substitutions can be seen from the 
structure of 6; (Eq.(8)). As screening length is different in electric and magnetic sectors the substitutions therefore 
involve different coefficients of mp and T for the transverse and the longitudinal case fia}. The longitudinal term 
after simplification like transverse one becomes, 


soft soft em? a? . ae . K 4 
if -aa = SBE rO Lile) — Lia(-e*)] + Ola’). (25) 








Again for the leading term we can write it in the following form, 
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where, 


w= 
— 
N 
ear 
SS 


a (=) = [I'(4) (Lis(—e7°) — Liz (—e®))]* 


The final expression for drag-coefficient then becomes, 
-mfi (T, ((B-w\\?_ 3x27, (Bn) 3 
"= E )aa\mp \ T T2773 \mp \ T 
6x 22/3 / T ‘i = 22 / T E 3 
ee ee Œ- T A Œ- ; (28) 
9711/3 \mp T 96E \mp T 


In the zero temperature limit the functions behave as h;(a) > |a| and g;(a@) > |a|. Hence, 7 in the extreme zero 
temperature limit becomes, 











o @ļE- u? 3x 2em} (zu) e?°|E — pl? (29) 


— 48TE T2n7/3E mp 96mpE 
This is the result for zero temperature plasma. Both the first and the second term here come from the transverse 
sector while the last piece emanates from the longitudinal interactions. The appearance of the second term with 
fractional power both in Eqs.(28) and clearly show that full contributions to 7 cannot be obtained by adding 
leading order contributions of the transverse and longitudinal photon exchange as the subleading terms of the former 
is larger than the leading order contribution of the later. This observation, in connection to the evaluation of Fermion 
self energy was first noted in and was overlooked in fal {16, . The zero temperature leading order contributions 
for | and t part are, however, consistent with our previous calculation reported in [18]. It is needless to mention here 
that such characteristic feature, also known as non-Fermi liquid behavior, can be attributed to the absence of the 
magnetostatic screening as noted in the introduction. 

In Fig.(2) we have plotted 7 versus energy of the incoming fermion in the small temperature (T/T; << 1) region 
where Ty = u/kpg is the Fermi temperature. From the figure it is evident that with increasing T/T, 7 decreases. 
This trend is consistent with what one finds for the fermionic damping rate at small temperature 

So far we have not discussed about the hard sector and tacitly assumed that the entire contribution to 7 in the 
relevant domains i.e. for small and zero temperature, come from the soft photon exchange. This, for degenerate 
plasma, is indeed so, as demonstrated explicitly in ig. In Eg, it was shown that the leading order of the hard sector 
fails to contribute to 7 at least up to O(e?). As, in the present work, on the other hand we go beyond the leading 
order, in principle, one should calculate the NLO part for the hard sector as well and see if the hard sector contributes 
to the drag and diffusion co-efficients in this case. But such an explicit calculation has not been done here. One can 
justify this omission on the ground that, for the soft sector, we see even after the inclusion of the NLO corrections 
no intermediate cut-off (q*) dependent term appear up to O(e?). Therefore, in the spirit of our previous work ig, 
we conclude up to this order the entire contribution comes from the soft sector providing indirect justification of this 
omission. The finite temperature NLO calculation can shed further light on this issue Bad. 
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FIG. 2: (Color online) The next to leading order drag coefficient at T/T; = 1/3 (dotted curve), T/T; = 1/2 (dashed curve), 
T/T; = 3/5 (dash dotted curve). 


B. Diffusion coefficient when |E — u| ~ T 


Along with 7, momentum diffusion coefficient (B;;), ig, p3, is another relevant quantity to study the 
equilibration of a fermion propagating in the plasma. For Coulomb plasma ņ and the longitudinal momentum diffusion 
coefficient (B) are related via Einstein’s Relation (ER). In this section we study the nature of longitudinal diffusion 
coefficient in the low temperature region. In the soft region the expression looks like, 





soft e2 q“ 1 
E) ~ an h daa’ f dzz? (1 +n(qx) — A(E — u — qz)) 
x {p(qz,q) + (1 = x°)p:(qx,q)}. 
(30) 


First, we calculate the transverse photon contribution then the longitudinal one. For the transverse photon propa- 
gator we proceed along the same line of the previous subsection and find, 
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The above expression is written in the following form, 
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After the magnetic part we derive the expression of the electric part. In case of electric term one finds, 


e = Seb E N sen (34) 
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e (==) = [F(5)(Li(-e7°) — Li4(—e°))] 


Finally, we obtain the expression for longitudinal momentum diffusion-coefficient as, 


sedeli /T (OSD ar N 
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This expression is polylogarithmic in nature and also contains fractional power in |E — u|. This fractional power 


indicates the deviation from Fermi-liquid behavior. This departure can also be seen in the zero temperature case. 
Hence, the final expression for B in the extreme zero temperature limit becomes, 
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ai 
= Saw orb (Ea)? EE A a, pi 


720 3377/3 mp 128mp 
The first two terms in the last equation correspond to the transverse contribution and the remaining third term comes 
from the longitudinal interaction. The expression for longitudinal diffusion coefficient has been already obtained in 
ag. Like 7, in B also we find that the subleading transverse part is greater than the leading longitudinal contribution. 
We see from the Fig.(8) that nature of the curve for the diffusion coefficient is same as that of 7 as shown in the 
previous subsection. 


II. CONCLUSION 


In this paper we have calculated the fermionic drag and diffusion coefficients in a relativistic plasma both at zero 
and small temperature by retaining terms beyond the leading contributions. It is seen that the subleading terms of 
the transverse sector, which appear with fractional power, are larger than the leading terms coming from the exchange 
of soft longitudinal photons or in other words we show that the leading order contributions to the drag and diffusion 
coefficients in ultradegenerate plasma cannot be obtained just by adding the leading order contributions coming from 
each of these sectors. Both the appearance of the fractional power and dominance of the transverse sector are related 
to absence of the magnetostatic screening or the singular behavior of the fermion self-energy near the Fermi surface. 
Furthermore, we find that the contributions coming from the hard sectors are suppressed and the entire physics is 
dominated by the soft excitations. This is a clear departure from the finite temperature case where both the hard and 
the soft part contribute at the same order. As a last remark, we note that here the entire calculation has been done 
for QED plasma. It would be interesting to extend the present calculation for QCD matter, which might be tricky 
due to the existence of triple gluon vertex and possible magnetic screening in the QCD sector. 
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FIG. 3: (Color online) The next to leading order diffusion coefficient at T/T; = 1/3 (dotted curve), T/T; = 1/2 (dashed curve), 
T/T; = 3/5 (dash dotted curve). 
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